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Abstract 



'sj" ' The existence of the three most massive clusters of galaxies observed so far at z > 0.5 is 

^ ■ used to constrain the mass density parameter of the universe, fl, and the amplitude of mass 

^ _ fluctuations, cxg. We find ri=0.2loj, and (Ts= 1.2lo;| (95 %). We show that the existence 

' of even the single most distant cluster at z=0.83, MS1054-03, with its large gravitational 

^ I lensing mass, high temperature, and large velocity dispersion, is sufficient to establish pow- 

' erful constraints. High-density, Q=l {as— 0.5-0.6) Gaussian models are ruled out by these 

^ ' data 10^^ probability); the Q=l models predict only ~ 10~^ massive clusters at z > 0.65 

^ ■ (~ lO^'^ at 2 > 0.5) instead of the 1 (3) clusters observed. 

subject headings .'cosmology : observation cosmology : theory galaxies : cluster : 



S ■ general galaxies : evolution large scale structure of universe 



1. Introduction 

The observed present-day abundance of rich clusters of galaxies places a strong con- 
straint on cosmology: crgfi*^'^ ~ 0.5, where as is the rms mass fluctuations on 8 h~^ Mpc 
scale, and Q is the present cosmological density parameter (Henry & Arnaud 1991, Bahcall 
& Cen 1992, White et al. 1993, Eke et al. 1996, Viana & Liddle 1996, Pen 1997, Kitayama 
& Suto 1997). This constraint is degenerate in Vt and ag; models with ^2=1, as^O.b are 
indistinguishable from models with f2~0.25, (Xg^l. (A as— 1 universe is unbiased, with mass 
following light on large scales; o"8~0.5 implies a mass distribution wider than light). 

The evolution of cluster abundance with redshift, especially for massive clusters, breaks 
the degeneracy between VL and as (see, e.g., Peebles et al. 1989, Oukbir & Blanchard 1992, 
1997, Eke et al. 1996, Viana & Liddle 1996, Carlberg et al. 1997, Bahcall et al. 1997, Fan 
et al. 1997, Henry 1997). The evolution of high mass clusters is strong in =1, low-crg 
(biased) Gaussian models, where only a very low cluster abundance is expected at z >0.5. 
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Conversely, the evolution rate in low-Jl , high-ug models is mild and the cluster abundance 
at z >0.5 is much higher than in 0=1 models. In Bahcall et al. (1997) and Fan et al. (1997) 
we used the CNOC cluster sample (Carlberg et al. 1997a,b, Luppino & Gioia 1995) to z <^ 
0.5 - 0.8 (with measured masses to ^ 0.5) to decouple and erg: we found ^2^0. 3 ±0.1 
and (Tg~0.83 ±0.15, consistent with Carlberg et a/.(1997a). The evolution rate, and the 
distinction among cosmological models, increases with cluster mass and with redshift: in Vt 
—1, low-CTg models, very massive clusters are not expected to exist at high redshifts. 

In the present paper wc extend the prcvioTis studies to larger mass and higher redshift 
clusters, using the three most massive clusters observed so far at high redshifts {z ~0.5-0.9) 
to independently constrain Vt and Ug. The clusters discussed in this paper are the three 
most massive distant clusters from the EMSS/CNOC sample used above, with masses larger 
by a factor of ~ 2 than the mass-threshold used previously (Evrard 1989, Bahcall et al. 
1997, Fan et al. 1997, Carlberg et al. 1997a). Rehably measured masses are now available 
for these clusters from gravitational lensing, temperatures, and velocity dispersions, not 
previously available in the above studies. Strong Sunyaev-Zel'dovich decrements have also 
been observed for these clusters, further suggesting that these are massive clusters with large 
amount of gas. The three clusters have the highest masses (from weak lensing observations), 
the highest velocity dispersions {a^ ;^1200 km <s~^), and the highest temperatures (T;^8 kev) 
in the z > 0.5 EMSS survey (§2). Therefore, they provide a strong constraint on cosmology. 
We discuss the cluster data in §2 and the cosmological implications in §3. A Hubble constant 
of Ho = 100 h km s-^Mpc"^ is used. 



2. Massive Distant Clusters 

We use the three most massive luminous clusters from the six 2; > 0.5 X-ray clusters in 
the EMSS Survey (Henry et al. 1992, Luppio & Gioia 1995). Extensive data are now available 
for these three clusters : temperatures (T 8 kev), velocity dispersions {a^ ^ 1200 km s~^), 
and masses (M(< 0.5h~-'^Mpc) 5 x lO-'^^h^^M©, from weak gravitational lensing observed 
for two of the three clusters). A summary of the data is presented in Table 1. 

The mass of the most distant cluster, MS1054-03 z — 0.83, has been determined 
from weak lensing (Luppino & Kaiser 1997): Mgl(< 0.82h-^Mpc) ~ lO^^h^^M© (within 
R = 0.82h~^Mpc = 1.5h~^ Mpc-comovm^^radius). This mass assumes the background galax- 
ies are located at the most realistic redshift range of 2; ~ 1.5—3 (Luppino & Kaiser 1997), and 
takes account of the surface mass zero point (Luppino & Kaiser 1997) and the de-projection 
to spherical mass within R (as given by Hjorth et al. 1997). The uncertainty in this mass 
estimate is ±30%. (If the background galaxies are located at a mean redshift of z ^ 1.5, the 
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cluster mass is even larger, by 35 %; Luppino & Kaiser 1997.) The evolution analysis pre- 
sented in this paper requires only a best-estimate of a THRESHOLD MASS for the clusters 
(within a given radius); i.e., cluster masses need only be > M(threshold) ± AM(threshold). 
As long as the clusters are above the given mass threshold, accurate knowledge of the cluster 
mass is unimportant. We determine a conservative best-estimate mass threshold later on. 

MS1054-03 is at least twice as massive as the Coma cluster (at the same radius). The 
independent mass estimates from the cluster temperature (12.3 kcv) and velocity dispersion 
(ar — 1360 km s~^, Donahue et al. 1997), assuming hydrostatic equilibrium, agree well 
with the Icnsing mass (Table 1). These estimates use the observed mean M(T) and M(o"r) 
mass relations (see e.g., Bahcall & Fan 1998; also Hjorth et al. 1997, Carlberg et al. 1997a,b): 
M(< Ih-^Mpc) ~ Tkev X lO^^h-^M© and (M < 1.5h-^Mpc) ~ 2.4^2 x 1.5h-i(Mpc)/G. Mass 
uncertainties are estimated from the observed uncertainties in T and Ur and the ~ 20% rms 
scatter in the observed relations. Cluster masses can be extended to other radii R using 
the observed cluster profile: M(< R) oc R for R ~ 0.5 - Ih-^Mpc and M(< R) cx R^-^^ for 
R ~ 1 — 2h^^Mpc (Carlberg et al. 1997b; also Fischer & Tyson 1996). The results are only 
weakly dependent on the exact profile slope; changing the slope by ±0.3 changes the results 
by ^ 15%. 

The other two clusters, MS0016-I-16 and MS0451-03, are the two most luminous X-ray 
clusters in the z > 0.5 EMSS sample; they have high velocity dispersions (a^ ^ 1200 km s~^) 
and high temperatures 8 kev). The mass of MS0016+16 has been determined from weak 
lensing (for R ^ 0.3h~^ Mpc; Small et al. 1995), and is consistent (when extrapolated to the 
same radius) with the mass obtained from the cluster velocity dispersion and temperature 
(Table 1). Strong Sunyaev-Zel'dovich decrements have also been observed for all three clus- 
ters (Carlstrom 1997). The fact that the independent mass determinations (from lensing, 
temperature, and velocity dispersion) for the three clusters all yield cluster masses that are 
consistent with each other (Table 1) clearly indicates the reliabihty of these mass estimates. 
The strong Sunyaev-Zel'dovich decrements detected in these clusters further support this 
conclusion, indicating clusters with large amounts of gas. The existence of substructure in 
any of these clusters will not significantly affect the results, as shown by the agreement of 
all three independent mass methods. Furthermore, direct simulations show a <^ 15% effect 
on the mass determination due to substructure in comparable overdensity clusters (Evrard 
et al. 1996). 

In the following cosmological analysis, only a mass threshold is needed for the clusters; 
the exact cluster mass above the threshold is unimportant. Based on the data presented in 
Table 1, the best estimated mean cluster masses are 10 ±2, 8.7 ±1.3, and ll±1.4xlO^^h~^M0 
for the three clusters (within 1.5h~^comoving Mpc). The uncertainties represent the formal 
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errors for the weighted average cluster mass as determined from the three independent meth- 
ods: lensing, temperature, and velocity dispersion. We adopt a conservative mass threshold 
of Mi.5_eom ^ X lO^^-^Mo (withiu 1.5h-^Mpc comoving radius). This threshold is 
safely near the low end of the observed cluster masses; the mass threshold uncertainties are 
included in the final results. If the cluster masses are higher, the derived constraints become 
even stronger. Furthermore, even if one of the clusters is eliminated from the sample, (e.g., 
if one, or even two cluster masses are below the threshold), the main conclusions remain 
essentially unchanged. 

The number density of massive (> 8x lO^'^^h^^MQ) clusters at 2; > 0.5 is determined using 
the T,{l/Vmax) method (e.g., Luppino & Gioia 1995, Donahue et al. 1997). We find n{z = 
0.5-0.65) ~ Sin X 10-^h=^Mpc-^ (2 clusters) and n(0.65-0.9) ~ 3.4lii x lO-^h^^Mpc"'^ (1 
cluster) for Q =1. The error bars represent 68 % confidence level assuming Poisson statistics 
and equal likehhood for each logn. For Q~0.2, n ~ 2 x 10~^h^Mpc~^ for each of the two 
bins. (The density may be larger by a factor of ^ 2 if any of the other z >0.5 EMSS clusters 
are as massive as these). How do these densities compare with the abundance of z ~ 
clusters? The abundance of Mi 5 > 8 x lO-'^'^h^^MQ clusters at 2; ~ is taken from several 
sources: the mass function of Bahcall & Cen (1993); the temperature function of Henry 
(1997) (see also Eke et al. 1996, Viana & Liddle 1996 and Pen 1997); and the temperature 
function of Edge et al. (1990). The temperature threshold corresponding to the above mass 
(within 1.5h~^comoving Mpc, at 2; ~ 0.05) is T ;^ 6.4 kev (see M(T,R) relation above). The 
local cluster abundances are all in the range ~ 1.3 — 2 x 10~^h^Mpc~^(xlO='=°-^). The velocity 
function of the EN ACS survey (Fadda et al. 1996, Mazure et al. 1996, Borgani et al. 1997), 
at the relevant Mi. 5 threshold (corresponding to ^ 1010 km s^"^) is ~ 3 X lO-VMpc"^ 
(based on 2 clusters in the complete z < 0.08 sample). 

The abundance of high mass clusters evolves slowly from 2;~0.8tO2;~0; the abundance 
at 2; ~ 0.8 is only a factor of ~ 5-10 times lower than at present. (The evolution may in fact 
be even slower since we conservatively included, as a lower limit, only the three best studied 
most luminous clusters). As we show below, the expected evolution rate in fl —1 (low as) 
models is ~ 10^ times faster. 



3. Constraining Q. and ug 

We compare the data with the expected evolution of cluster abundance using the Press- 
Schechter (1974; P-S) formalism, which describes the evolution of the abundance of bound 
objects that grow from random-phase Gaussian initial fluctuations. The P-S method yields 
results that are in good agreement with simulations (down to the simulation limit of ~ 
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10-%^Mpc-^;e.g., Eke et al. 1996, Fan et al. 1997, Pen 1997). P-S describes the evolution 
of the cluster mass function. 



where a{z^ M) is the linear theory rms mass density fluctuation in spheres of mass M at 
redshift z, 5c — 1.68 is the critical density contrast needed for collapse (weakly dependent 
on Q, Eke et al. 1996, Kitayama & Suto 1996), and p is the mean cosmic density. The mass 
refers to the virial mass of the system. The present rms mass fluctuation within a sphere of 
mass M, cto(M), relates to ag as cto(M) = a^^M/M^)-'^ - crgM-'^l]", where a = (n+3)/6, 
n is the slope of the power spectrum at ~ 8h~^Mpc, and Mg oc Q is the mean mass within a 
sphere of radius 8h~^Mpc. We use n ~ —1.4, as observed (corresponding to Oh ~ 0.25 for a 
CDM spectrum); the results are insensitive to n and a (changing Qhy < 20% for reasonable 
changes in n; Fan et al. 1997). 

We use three different (but not independent) approaches. First, we calculate the evo- 
lution of the abundance of Mi.5_com > 8 x lO^^h^^M© clusters by integrating eq. (1). We 
use the O-cts constraint placed by the present-day cluster abundance: crgO^-^^ = 0.53 ± 0.05 
(for A = 0) and (jgO^-^^ = 0.53 ± 0.05 (for Q + A = 1) (Pen 1997; similar results by Eke et 
al. 1996; see also Bahcall & Cen 1992, White et al. 1993, Viana & Liddle 1996, Kitayama 
& Suto 1997); this ensures that the proper normalization cig is used for any O. The virial 
mass in eq. (1) is converted to the observed Mi.5_com mass as follows (see Fan et al. 1997): 
the virial overdensity is numerically calculated for each O (ag) model (e.g. Eke et al. 1996, 
Oukbir & Blanchard 1997), thus yielding the virial radius for a given virial mass. The virial 
mass (i.e., the mass within the virial radius) is then scaled to the 1.5 h~^ Mpc comoving 
radius using the observed M(< R) oc R°-^^ profile (see §2). The results are insensitive to this 
transformation; changing the profile slope by ±0.3 changes the constraints on O by ^ 15%. 

Figure 1 shows the expected evolution of the cluster abundance as a function of redshift 
for different O, erg combinations that satisfy the present day cluster abundance. The model 
curves range from Q — 0.1 (dg ~ 1.7) at the top of the figure (fiattest, nearly no evolution) to 
O = 1 ((Tg ~ 0.5) at the bottom (steepest, strongest evolution). The difference between low 
and high Q models (i.e., high and low cig) is dramatic for these high mass clusters; the cluster 
abundance in O = 1 models is ~ 10^^ of the abundance in low Q models. The data exhibit 
only a slow, relatively flat evolution; this is expected only in low-Q models. The observed 
clusters, even just the single cluster at 2; ~ 0.83, represent ~ 10^ more clusters than expected 
in any O =1 (ag ~ 0.5-0.6) models such as Cold or Mixed Dark Matter (CDM or MDM). 
The results are similar for A =0 and A +Q —1 models (the latter show shghtly stronger 
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evolution; see Bahcall et al. 1997 and Fig. 2). 

Figure 2 shows the allowed range (68 % and 95 %) of Q permitted by the existence of 
the massive clusters aX, z — 0.5 — 0.9 (using the crg-Jl normalization above). The expected 
number of clusters with Mi.5_com > 8 x lO^^h^^M© in the two redshift bins in the current 
sample is presented as a function of (or cts), and compared with the two and one detected 
clusters, respectively. The existence of the two clusters ai z = 0.5 — 0.65 is only possible if 
Vt ^ 0.25to;J^, (Tg — 1-2 ± 0.3 (95%). As Q increases (cxs decreases), the expected number 
of clusters drops exponentially. For Q = 1 (ag ~ 0.5), only ~ 10"'^ clusters are expected 
compared with the 2 observed, representing a probability of ~ 10~^. Similar (but indepen- 
dent) results are obtained for the one cluster observed in the z = 0.65-0.9 bin (Fig. 2); 
the probabihty of finding one such cluster if Q = 1 is ~ 10~^. These independent results 
reinforce the above conclusions, ruling out models as very low probability (^ 10~^). 
This low probability (^ 10~^) applies for any mass-threshold M 7 x lO^^h^^Mg; the larger 
the mass, the lower the probability. The difference between A = and A +D, — 1 models is 
small (Fig. 2). 

A second method of analyzing the results, independent of the ag-ft normalization, was 
discussed by Fan et al. (1997). We showed that the evolution rate of cluster abundance 
depends exponentially on erg; if we ignore the normalization crsfi"'^ ~ 0.5 that is needed to 
match the "absolute" cluster abundance at z ~ and study only the rate of the evolution 
profile with z (i.e., its "flatness"), we can estimate directly, nearly independent of Q (as 
well as nearly independent of the exact shape of the power spectrum and Hq). Figure 3 
shows the cluster abundance ratio , n{z ~ 0)/n{z ~ 0.8), as a function of erg (for all fl 's), 
determined from P-S (eq. 1). The very strong evolution for low -erg is clearly distinguished 
from the nearly no evolution expected for Ug ~ 1 models (for any ). The slow observed 
evolution rate (Fig.3) requires a high erg value : erg ^ 0.9, independent of Q . When com- 
bined with the proper normalization of the present-day cluster abundance, a low ft value is 
obtained: Q ^ 0.35. Similar results are found for A = 0. 

A third method of analyzing the results is to determine an independent Q -erg relation 
at high redshift (using P-S) based entirely on the three high redshift clusters; this method 
is independent of the z=0 cluster abTindance. Figure 4 shows the allowed fl -cxg range 
determined by the high-z cluster abundance. The two high redshift bins (shown by the solid 
and dashed curves in Fig. 4) are consistent with each other. They satisfy 

(7gQ°-^^ ~ 0.8±0.1 (A = 0), (7gQ°-^^ ~ 0.8±0.1 (A + Q = l) [for z ~ 0.5 - 0.9] (2) 

The dependence on Q is weak. These relations imply a high erg normalization, erg ^ 0.8, 
for any ft < 1. The independent constraint provided by the present-day cluster abundance 
(§3) is also presented in Fig. 4. These two Q -erg relations overlap only for low-Jl , high-erg 
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values. The allowed parameter ranges (68 % and 95 % confidence level contours) are shown 

in Fig 4b; the probability distributions are derived from the likelihood of each Q — ag, model 
assuming Poisson statistics. We find = 0.22lo;j7 and (T8=l-18lo:22 (for A + Q = 1), and 
n=0.17t'^o^ot , as=l.l7t'^oil (for A = 0) (68 % limits). The 95% limits arc n = 0.22l|]:f , 
as = 1.18l[J:^2 (1^ + A = 1) and = 0.17l{]:f , = l-Ut^oH = 0). The resuhs are 
consistent (at ~ 1 a level) with those obtained previously by Bahcall et al. (1997) and Fan 
et al. (1997) using lower mass clusters, and with the low-Jl conclusion of Luppino & Kaiser 
(1997), Henry (1997), and Donahue et a;.(1997). 

In summary, wc find Cl = 0.21q;'| and = 1.2^q'^ (95% C.L.). The high cg value suggests 
a nearly unbiased universe, where mass approximately follows light on large scales. The error- 
bars include an estimated uncertainty in the cluster mass threshold ranging from Mi.5_com — 
7 to 10 X lO^^h^^M©. The existence of even the single massive cluster at 2;=0.83 places 
similar stringent constraints. Using a lower mass threshold of Mi.5_com = 7 x lO^^h^^MQ 
increases by 30 %. The results are insensitive to the power spectrum and to Hq; fixing h 
= 0.65 ± 0.15 (corresponding to n ~ — 1 for Q = 1, instead of n ^ —1.4) lowers as by 10% 
and increases Q by 20%. To be consistent with O =1 models (such as SCDM or MDM, as — 
0.5-0.6) all three cluster masses need to be unacceptably low. Mi 5_com ~ 2 — 3 x 10"'^^h~^Mpc 
(i.e. comparable to poor clusters of T ~ 2 - 3 kev, instead of the observed T ~ 8 — 12 kev); 
this is excluded by the lensing, velocity dispersion, temperature, and Sunyaev-Zel'dovich 
data. 

We thank J. Hughes, G. Luppino, J. P. Ostrikcr, J. P. E. Peebles, D. N. Spcrgcl and M. 
A. Strauss for helpful discussions. The work is supported by NSF grant AST93- 15368. 
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Table 1. Cluster Properties 



Cluster 


z 


ivr(a) 
-'^0.5 


kev 




-f^l.5 — com 




,)(10i^h- 












(km s-i) 


(h-^Mpc) 


Mg.l. M(T) 


M(C7) 


< M > 


MS1054-03 


0.83 


50 ± 10 




1360 ± 450 


0.82 


^10 ±3 10 ±3 


10 ±6 


10 ±2 


MS0016+16 


0.55 


66 ± 16 


8 ± 1 


1234 ± 128 


0.97 


9±4 8±2 


9±2 


8.7 ±1.3 


MS0451-03 


0.54 


47 ± 5 


10.4 ± 1.2 


1371 ± 105 


0.97 


10 ±2 


12 ±2 


11 ±1.4 



a) Central galaxy density count (Luppino & Gioia 1995). {No.5{Coma) = 28; Bahcall 1977) 

b) Cluster temperature (Mushotzky & Scharf 1997, Donahue 1996, Donahue et al. 1997; 90 % error bars) 

c) Cluster velocity dispersion (Carlberg et al. 1997b for z ~ 0.5 clusters, 68 % error bars; Donahue et al. 
1997 for z ~ 0.8 cluster, 90 % error bar) 

d) Cluster mass (within Rcom = 1.5h~-'^Mpc, column 6) as determined from gravitational lensing (Mg.l.; 
Small et al. 1995, Luppino & Kaiser 1997, Hjorth et al. 1997), and from the cluster temperature and velocity 
dispersion of column 4-5 (M(T), M(cr); §2). The weighted mean estimated mass for each cluster based on 
these data, < M >, is also listed. Only a threshold mass is needed for the evolution analysis; we adopt, 
conservatively, Mi.5_coin > ^ti x lO^^h~^M0 for these clusters (§2). 
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Figure 1. Observed vs. model evolution of cluster abundance for M(< Rcom = 1.51i~^Mpc) > 

8 X lO^^h"^M0 clusters. The P-S model expectations (curves) are shown as a function of 
^(ag) (for O + A = l;A = models yield similar results). The z ~ 0.6 and 0.8 data points 
represent 2 and 1 clusters, respectively (the most massive EMSS clusters; §2); filled circles 
are for 0=1, open circles 9. ^ 0.2 (A = 0), triangles ^ ~ 0.2 (A = 0.8). The 68% and 95% 
error bars are shown. (For 2; ~ 0, see §2). 
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Figure 2. Observed vs. expected (P-S) number of clusters in the survey in each of the 
two z > 0.5 bins (Fig 1), presented as a function of Q (or erg ). Sohd hues represent the 
z = 0.65 — 0.90 bin (1 observed cluster); dashed lines represent the z = 0.5 — 0.65 bin (2 
observed clusters). The 68 % and 95 % confidence ranges are marked by the dark and open 
regions near Q ~ 0.2, Ug ~ 1. The probability oi Q = 1 models being consistent with either 
data bin is ~ 10~^. 
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Figure 3. Cluster abundance ratio, n{z = 0.8)/n{z = 0), vs. erg from Press-Schechter (solid 
curve, for mean of all Q's) for Mi.5_com > 8 x lO^'^M© clusters. Filled circles represent Q's 
from 0.2 to 1 (bottom to top). (Dashed curves represent the mass threshold range of 7 to 
10 xIO^^Mq, top and bottom, respectively). The data (Fig 1, §2) are shown by the shaded 
region (68 % level). Similar results are obtained for A = (with weaker dependence on Q). 




^8 

Figure 4. Constrainting Q and as- The z = band represents the present-day cluster 
abundance constraint (crgQ^-^ ~ 0.5; Pen 1997 dot-dashed curve, Eke et al. 1996 dotted 
curve). The z — 0.6, 0.8 bands (dashed and sohd curves) are the constraints provided by the 
present high redshift clusters (crgfi^'^^ ~ 0.8; eq 2.) (68%). The allowed parameter ranges 
(68 % and 95 %) are shown (bottom figure). Similar results are obtained for A = (§3). 
The COBE normahzation (Bunn & White, 1997) are also shown (open and A models). 



